{smcl}
{com}{sf}{ul off}{txt}{.-}
       log:  {res}descriptive statistics.smcl
  {txt}log type:  {res}smcl
 {txt}opened on:  {res} 9 Jan 2011, 17:24:56

{com}. do "C:\Users\Piotr\AppData\Local\Temp\STD06000000.tmp"
{txt}
{com}. **** DESCRIPTIVE STATISTICS ON MIPT with "domestic" Iraq
. 
. clear
{txt}
{com}. set memory 20m
{txt}(20480k)

{com}. 
. ************* "ALL" **********
. 
. use MIPT_IraqAll.dta, clear
{txt}
{com}. 
. * All attacks in all income classes by quarter over 40 yrs
. twoway line mAll Quarter, title("Panel 1: All incidents") legend(off) ttick(1968q1(4)2007q1) tlabel(1968q1(12)2007q1, format(%tqCCYY) labsize(3)) tmtick(1968q3(4)2007q4) ytitle("Incidents per quarter") ylabel(0(100)1100, angle(horizontal) labsize(3)) ytick(0(50)1100) tline(177) saving(IraqAll_all.gph, replace)
{res}{txt}(file IraqAll_all.gph saved)

{com}. 
. * Attacks with fatalities in all income classes by quarter:
. twoway line mFatal Quarter, ttick(1968q1(4)2007q1) tlabel(1968q1(12)2007q1, format(%tqCCYY) labsize(3)) tmtick(1968q3(4)2007q4) ytitle("Incidents per quarter") ylabel(0(100)800, angle(horizontal) labsize(3)) saving(IraqFatal_all.gph, replace)
{res}{txt}(file IraqFatal_all.gph saved)

{com}. 
. * Attacks with injuries in all income classes by quarter:
. twoway line mInj Quarter, ttick(1968q1(4)2007q1) tlabel(1968q1(12)2007q1, format(%tqCCYY) labsize(3)) tmtick(1968q3(4)2007q4) ytitle("Incidents per quarter") ylabel(0(100)500, angle(horizontal) labsize(3)) saving(IraqInj_all.gph, replace)
{res}{txt}(file IraqInj_all.gph saved)

{com}. 
. * Attacks with any casualties (fatalities and injuries) in all income classes by quarter:
. twoway line mCasualty Quarter, title("Panel 3: All casualty incidents") legend(off) ttick(1968q1(4)2007q1) tlabel(1968q1(12)2007q1, format(%tqCCYY) labsize(3)) tmtick(1968q3(4)2007q4) ytitle("Incidents per quarter") ylabel(0(100)1100, angle(horizontal) labsize(3)) ytick(0(50)1100) tline(179) saving(IraqCasualty_all.gph, replace)
{res}{txt}(file IraqCasualty_all.gph saved)

{com}. 
. * Attacks against US targets:
. twoway line mUS_Target Quarter, title("Panel 4: All incidents against a U.S. Target") legend(off) ttick(1968q1(4)2007q1) tlabel(1968q1(12)2007q1, format(%tqCCYY) labsize(3)) tmtick(1968q3(4)2007q4) ytitle("Incidents per quarter") ylabel(0(100)1100, angle(horizontal) labsize(3)) ytick(0(50)1100) tline(177) saving(IraqUStarget_all.gph, replace)
{res}{txt}(file IraqUStarget_all.gph saved)

{com}. 
. 
. drop mInj mFatal mUSFatal mUSInj mUSCasualty
{txt}
{com}. tsset Quarter
{res}{txt}{col 9}time variable:  {res}{col 25}Quarter, 1968:1 to 2007:4
{txt}{col 17}delta:  {res}1 quarter
{txt}
{com}. save IraqAll_all.dta, replace
{txt}file IraqAll_all.dta saved

{com}. clear
{txt}
{com}. 
. 
. ************** "MICs" *****************
. *Adding "domestic" incidents in Iraq changes numbers only for MICs - no need to present LICs and HICs
. 
. clear
{txt}
{com}. use MIPT_IraqMIC.dta, clear
{txt}
{com}. 
. * All attacks in MICs by quarter over 40 yrs
. twoway line mAll Quarter, title("Panel 2: MIC incidents") legend(off) ttick(1968q1(4)2007q1) tlabel(1968q1(12)2007q1, format(%tqCCYY) labsize(3)) tmtick(1968q3(4)2007q4) ytitle("Incidents per quarter") ylabel(0(100)1100, angle(horizontal) labsize(3)) ytick(0(50)1100) tline(177) saving(IraqAll_MIC.gph, replace)
{res}{txt}(file IraqAll_MIC.gph saved)

{com}. 
. * Attacks with fatalities in MICs by quarter:
. twoway line mFatal Quarter, ttick(1968q1(4)2007q1) tlabel(1968q1(12)2007q1, format(%tqCCYY) labsize(3)) tmtick(1968q3(4)2007q4) ytitle("Incidents per quarter") ylabel(0(50)800, angle(horizontal) labsize(3)) saving(IraqFatal_MIC.gph, replace)
{res}{txt}(file IraqFatal_MIC.gph saved)

{com}. 
. * Attacks with injuries in MICs by quarter:
. twoway line mInj Quarter, ttick(1968q1(4)2007q1) tlabel(1968q1(12)2007q1, format(%tqCCYY) labsize(3)) tmtick(1968q3(4)2007q4) ytitle("Incidents per quarter") ylabel(0(50)500, angle(horizontal) labsize(3)) saving(IraqInj_MIC.gph, replace)
{res}{txt}(file IraqInj_MIC.gph saved)

{com}. 
. * Attacks with any casualties (fatalities and injuries) in MICs by quarter:
. twoway line mCasualty Quarter, title("Panel 3: MIC casualty incidents") legend(off) ttick(1968q1(4)2007q1) tlabel(1968q1(12)2007q1, format(%tqCCYY) labsize(3)) tmtick(1968q3(4)2007q4) ytitle("Incidents per quarter") ylabel(0(100)900, angle(horizontal) labsize(3)) saving(IraqCasualty_MIC.gph, replace)
{res}{txt}(file IraqCasualty_MIC.gph saved)

{com}. 
. * Attacks against US targets in MICs:
. twoway line mUS_Target Quarter, title("Panel 3: MIC incidents against a U.S. target") legend(off) ttick(1968q1(4)2007q1) tlabel(1968q1(12)2007q1, format(%tqCCYY) labsize(3)) tmtick(1968q3(4)2007q4) ytitle("Incidents per quarter") ylabel(0(10)100, angle(horizontal) labsize(3)) saving(IraqUStarget_MIC.gph, replace)
{res}{txt}(file IraqUStarget_MIC.gph saved)

{com}. 
. drop mInj mFatal mUSFatal mUSInj mUSCasualty
{txt}
{com}. tsset Quarter
{res}{txt}{col 9}time variable:  {res}{col 25}Quarter, 1968:1 to 2007:4
{txt}{col 17}delta:  {res}1 quarter
{txt}
{com}. save IraqMIC_all.dta, replace
{txt}file IraqMIC_all.dta saved

{com}. 
. clear
{txt}
{com}. 
. 
. ************* "ONLY IRAQ (intl.+domestic after invasion)" **********
. 
. use MIPT_only Iraq ALL for regression.dta, clear
{txt}
{com}. 
. * All attacks in Iraq by quarter over 40 yrs
. twoway line mAll Quarter, title("Panel 2: All Iraq incidents (extended series)") legend(off) ttick(1968q1(4)2007q1) tlabel(1968q1(12)2007q1, format(%tqCCYY) labsize(3)) tmtick(1968q3(4)2007q4) ytitle("Incidents per quarter") ylabel(0(100)1100, angle(horizontal) labsize(3)) ytick(0(50)1100) tline(177) saving(IraqONLY_all.gph, replace)
{res}{txt}(file IraqONLY_all.gph saved)

{com}. 
. * Attacks with any casualties (fatalities and injuries) in all income classes by quarter:
. twoway line mCasualty Quarter, title("All Iraq casualty incidents") legend(off) ttick(1968q1(4)2007q1) tlabel(1968q1(12)2007q1, format(%tqCCYY) labsize(3)) tmtick(1968q3(4)2007q4) ytitle("Incidents per quarter") ylabel(0(100)1100, angle(horizontal) labsize(3)) ytick(0(50)1100) tline(177) saving(Only IraqCasualty_all.gph, replace)
{res}{txt}(file Only IraqCasualty_all.gph saved)

{com}. 
. 
. ************* "MICs excluding Iraq" **********
. 
. use MIPT_MICs withoutIraq for regression.dta, clear
{txt}
{com}. 
. * All attacks in MICs excluding Iraq by quarter over 40 yrs
. twoway line mAll Quarter, title("Panel 4: All MICs incidents excluding Iraq") legend(off) ttick(1968q1(4)2007q1) tlabel(1968q1(12)2007q1, format(%tqCCYY) labsize(3)) tmtick(1968q3(4)2007q4) ytitle("Incidents per quarter") ylabel(0(20)150, angle(horizontal) labsize(3)) ytick(0(10)150) saving(MICs_excl_Iraq_all.gph, replace)
{res}{txt}(file MICs_excl_Iraq_all.gph saved)

{com}. 
{txt}end of do-file

{com}. clear

. do "C:\Users\Piotr\AppData\Local\Temp\STD06000000.tmp"
{txt}
{com}. *** TABLE WITH DESCRIPTIVE / SUMMARY STATISTICS ***
. *with domestic Iraq
. clear
{txt}
{com}. set memory 20m
{txt}(20480k)

{com}. 
. *** ALL ***
. 
. use MIPT_IraqAll.dta
{txt}
{com}. 
. *All incidents*
. 
. *2003:3 2007:4
. sum mAll if(Quarter >173)

{txt}    Variable {c |}       Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 56}
        mAll {c |}{res}        18    569.9444    339.9509         82       1120
{txt}
{com}. 
. *Casualty incidents*
. 
. *2003:3 2007:4
. sum mCasualty if(Quarter >173)

{txt}    Variable {c |}       Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 56}
   mCasualty {c |}{res}        18    442.7778    287.7966         39        961
{txt}
{com}. 
. *Incidents against a U.S. target*
. 
. *2003:3 2007:4
. sum mUS_Target if(Quarter >173)

{txt}    Variable {c |}       Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 56}
  mUS_Target {c |}{res}        18    540.3889    334.7642         56       1103
{txt}
{com}. 
. *Casualty incidents against a U.S. target*
. 
. *2003:3 2007:4
. sum mUSCasualty if(Quarter >173)

{txt}    Variable {c |}       Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 56}
 mUSCasualty {c |}{res}        18    432.1111    286.5589         36        950
{txt}
{com}. 
. clear
{txt}
{com}. 
. ************** "MICs" *****************
. 
. clear
{txt}
{com}. use MIPT_IraqMIC.dta
{txt}
{com}. 
. *All incidents*
. 
. *2003:3 2007:4
. sum mAll if(Quarter >173)

{txt}    Variable {c |}       Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 56}
        mAll {c |}{res}        18    545.8889    335.3365         56       1106
{txt}
{com}. 
. *Casualty incidents*
. 
. *2003:3 2007:4
. sum mCasualty if(Quarter >173)

{txt}    Variable {c |}       Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 56}
   mCasualty {c |}{res}        18    433.9444    287.0107         35        953
{txt}
{com}. 
. *Incidents against a U.S. target*
. 
. *2003:3 2007:4
. sum mUS_Target if(Quarter >173)

{txt}    Variable {c |}       Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 56}
  mUS_Target {c |}{res}        18       537.5     335.329         50       1100
{txt}
{com}. 
. *Casualty incidents against a U.S. target*
. 
. *2003:3 2007:4
. sum mUSCasualty if(Quarter >173)

{txt}    Variable {c |}       Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 56}
 mUSCasualty {c |}{res}        18    430.4444    286.5795         35        948
{txt}
{com}. 
{txt}end of do-file

{com}. log close
       {txt}log:  {res}descriptive statistics.smcl
  {txt}log type:  {res}smcl
 {txt}closed on:  {res} 9 Jan 2011, 17:30:49
{txt}{.-}
{smcl}
{txt}{sf}{ul off}